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Abstract. We study orientational order, subject to thermal fluctuations, on a 
fixed curved surface. We derive, in particular, the average density of zeros of 
Gaussian distributed vector fields on a closed Riemannian manifold. Results are 
compared with the density of disclination charges obtained from a Coulomb gas 
model. Our model describes the disordered state of two dimensional objects with 
orientational degrees of freedom, such as vector ordering in Langmuir monolayers 
and lipid bilayers above the hexatic to fluid transition. 



PACS numbers: 02.40.-k, 87.16.Dg, 61.30.Jf 

1. Introduction 

In several areas of statistical physics and condensed matter, a great deal of progress 
has been achieved by focusing on the physics of topological defects, ignoring other 
degrees of freedom. The Kosterlitz-Thouless transition, describing the destruction of 
orientational order in thin films, is a particularly important example. The transition 
is viewed as one where defect pairs unbind and proliferate, destroying the (quasi)- 
long-range order Q. Another crucial example is that of typc-II superconductors in 
a magnetic field, where the important physics is encoded in the properties of vortex 
lines |2[|. Even in the absence of a magnetic field, the formation and growth of vortex 
loops can be used to explain the form of the voltage vs. current relation ||. Yet 
a third example is found in the physics of orientational order in membranes, where 
many workers @, H have found it fruitful to focus on the properties of topological 
defects to understand the low-temperature physics. 

In all of these examples, topological defects are used primarily to understand the 
low-temperature behavior. For example, in the case of superconductors, the vortex 
line description is used primarily to understand the behavior below T c (or H C 2), rather 
than to understand the properties of the normal metal phase at higher temperatures. 
Indeed, topological defects are a more natural description at low temperatures, where 
it is very costly to excite the order parameter away from its average value. Although 
themselves energetically costly, topological defects are the minimal-energy way of 
satisfying a constraint of the system - e.g., the curvature of a membrane or the 
penetration of a magnetic field into a superconductor. 

At higher temperatures, the description becomes less natural. Order parameter 
fluctuations become much less costly, and hence the fluctuations observed in thermal 
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equilibrium become more violent. At sufficiently high temperatures, the broken 
symmetry associated with the low-temperature phase is restored, and the average 
value of the order parameter is zero. Above this temperature, the order parameter 
behaves approximately like a Gaussian random variable. In this case, it seems that the 
description in terms of topological defects will be insufficient to describe the physics, as 
it only encodes the locations (and signs) of the zeros of the order parameter, classifying 
as irrelevant any fluctuations that the order parameter undergoes between these zeros. 

Despite these objections, a description in terms of topological defects actually 
approximates the high temperature behavior of such systems quite well. For the case 
of thin films, Halperin H showed that the density of topological defects (i.e. zeros 
of the order parameter) that one obtains from a Gaussian order parameter are much 
the same as one obtains from a Coulomb gas model that allows only the topological 
defects as degrees of freedom. In the case of superconductors, Lehrer and Nelson || 
have shown that above H C 2, a Gaussian approximation to the Ginzburg-Landau free 
energy predicts approximately the same distribution of vortex loops and lines as does 
the London theory, which is purely a description in terms of topological defects. 

In this paper, we focus on the case of membranes. In particular, we examine the 
density of topological defects under the approximation that the free energy is Gaussian. 
This approximation will be valid at high temperatures. We compare these results to 
results obtained from a model which focuses only on the topological defects and their 
interactions, namely, a Debye-Huckel approximation to a Coulomb gas model. 

One way to probe the high-temperature properties of a topologically spherical 
surface is with light scattering experiments on lipid vesicles Q ■ In the case of lipid 
bilayers, the source of the orientational order parameter is the vector between a lipid 
head and a neighboring head, and the order parameter describes the tendency of 
the lipids to achieve hexatic order at low temperatures. (The hexatic order and the 
transition to the disordered phase were studied in free standing liquid-crystal films 
using light and X-ray scattering [pT[.) However, since this order parameter is invariant 
under rotations by 60° (because most lipids have six neighbors on average) , it is slightly 
different from the case we consider, where the order parameter is only invariant under 
rotation by integer multiples of 360°. Nevertheless, we expect much of what we derive 
here to apply to these systems. 

A system that is closer to what we consider here is that of tilted Langmuir 
monolayers [|f2"| , which consist of a monolayer of lipids or amphiphilcs on a liquid 
surface, e.g. the surface of a water droplet. The surfactant molecules have a tendency 
to orient themselves so that the hydrocarbon chain is at a tilt away from the normal 
to the surface. The projection of the direction of the polymer chain onto the surface 
forms an orientational order parameter which is exactly of the form that we consider 
in this paper, as illustrated in figure [|. A similar situation may be obtained in a 
lipid bilayer when the lipids tend to tilt away from the normal to the surface |l3|, [b| , 
providing yet another source for an orientational order parameter. 

The rest of the paper is organized as follows. In section ||, we write down 
a continuum model for the orientational order parameter that has the correct 
symmetries and should describe both the low- and the high-temperature physics 
of the orientational order parameter, namely, a 0(2)— Ginzburg-Landau theory in 
curved space. From this, we will derive a Coulomb gas model, and use it to calculate 
the density of disclination defects in the Debye-Hiickel approximation in the high- 
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temperature limit. We obtain 

2np = K+ * AK + 0(x- 2 ), (1.1) 

where p is the density of defects, K is the curvature of the surface, A is the Laplace- 
Beltrami operator, x is the fugacity of the topological defects, and Ka describes the 
interaction strength of the defects. 

In section |^, we present the main results of our paper. We approximate our 
model by neglecting the nonlinear terms, valid at high temperatures, and calculate 
the density of defects for an arbitrarily curved surface. For simplicity, we restrict our 
scope to closed surfaces - namely surfaces that are topologically equivalent to spheres, 
tori, etc. We especially focus on the case where the surface is topologically equivalent 
to a sphere, both for calculational ease and because we expect this class of (closed) 
surfaces to be the most easily amenable to experiments. Because the calculation is 
fairly technical, we first review Halperin's calculation of the density of topological 
defects that we expect to see in flat space as a "toy model" for the problem in curved 
space. We then proceed to the calculation in curved space, obtaining 

AK 

27t P = K + —— + 0(t- 2 ), (1.2) 
IzttZt 

where r measures the deviation from the critical temperature, Z^ 1 = 2n/ log(l/(a 2 T)), 
and a is a short distance cutoff. Instead of using a gauge-field representation of 
the orientational order parameter we deal with a manifestly gauge-invariant picture. 
Employing a special symmetry of the model we can express the director field through 
simple scalar fields and solve it exactly in a high-temperature expansion. 

The result for the defect density in curved space is equivalent to the Coulomb 
gas/Debye-Hiickel result above, provided we identify ttKax = 3Zt. This confirms 
the validity of the Coulomb model even for high temperatures at this level of 
approximation. Deviations do, however, begin to show up at 0(r -2 ), as we shall 
show below. 



2. Model and Debye-Hiickel theory 

We concentrate on the case of purely in-plane orientational order and therefore use a 
(two component) tangential vector field Ui(a) as the order parameter. To describe the 
physics of the surface, we rely on the language of Riemannian differential geometry, 
which ensures that results are independent of any particular coordinate system. (A 
concise introduction to differential geometry of surfaces can be found in reference pqj .) 
The order parameter lives on a closed two-dimensional Riemannian manifold with 
line element ds 2 = gijdo* da 3 , where g-y = gij(a 1 , a 2 ) is the metric tensor and 
a = (a 1 , a 2 ) are internal coordinates of the surface. Using this formalism, we write 
down a O (2)— invariant, statistical weight P[u] cx exp(— H/T) for the Uj— field, where 
H is the (mesoscopic) Hamiltonian and T is the temperature. The simplest such 
Hamiltonian H is the analogue of Ginzburg-Landau theory in flat space, namely 

^ = ^JdA (DiUjITv* + ruiu' + cfat**) 2 ) , (2.1) 

where dA = ^/g(<j) dodder 2 is the invariant area element, g is the determinant of 
the metric, Di is the covariant gradient, u % — g i0 Uj, g v = {g~ r )ij, c is the coupling 
constant, and the coupling constant for the gradient term is absorbed into the field 
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it, . Equation ( p.lj ) encodes the same physics as the free energy used by Park et al Q 
and by Evans p] for vector defects, although theirs appears in the gauge-field picture. 



The equivalence of the models is shown explicitly in Appendix C . 

For t below a certain critical r, this will be a "Mexican hat potential" ; however, 
we will be concerned primarily with the opposite case of high temperatures (t > 0). 
Before we specialize to this case, we look at some properties at smaller r. A critical 
temperature r c (the mean field value r c = gets renormalized due to fluctuations) 
separates the disordered state r > r c (high-temperature region) from the "ordered 
state" t < t c . We put "ordered state" in quotation marks, because a perfectly ordered 
state is impossible for certain manifold topologies. For example, on a sphere or any 
other surface with the same topology, a tangential vector field has at least two zeros 
(defects) Jl5| . This can be illustrated by attempting to comb a hedgehog or a hairy 
ball: there will be two places where the vector field is zero or has a singularity. 

To investigate this in more detail, we distinguish between two types of zeros. One 
type, called a "positive zero", is characterized by det(DiUj) > 0, while the other type, 
a "negative zero", has a saddle-like flow and is characterized by det(DiUj) < 0. See 
figure H for an illustration of these types of defects. Zeros with det (DiUj) = do 
not fall into this scheme; however, they will not show up in a statistical model as the 
probability to hit exactly det(DiUj) = vanishes. The number of positive zeros minus 
the number of negative zeros is a topological constraint and equal to 2(1 — 7), where 7 
is the number of handles of the (closed) surface, e.g. zero for a spherical topology and 
one for a torus We will show this theorem explicitly en route to our calculation. 
Thus, on surfaces other than tori, the low-temperature phase necessarily has defects, 
unlike in flat space, where the ground state is defect free. 

As in flat space, the properties of the low temperature phase are determined by low 
energy Goldstone modes ("spin waves"), which prevent true long ranged correlations. 
Instead, one finds an algebraic decay of the correlations (quasi- long-ranged order). 
Besides the spin waves, thermally excited defects persist. Integration over the spin 
waves results in a Coulomb gas model for these defects (zeros) M, where the defects 
carry a charge proportional to their index q = sign det (DiUj) and a core energy 
(chemical potential Q ) . The interaction energy of the defects reads 

f = ^ J dA J dA' (2irp - K) a G{a, o>) (2irp - K)„, , (2.2) 

where G(a, a') denotes the Green's function of the negative of the Laplace-Beltrami 
operator —A = —g^Didj | and Ka is the coarse-grained effective coupling between 
the defects and therefore depends on the temperature |], [l7j] . p is the defect density 
p(a) = qi 5 c {a, Ci), where ai are the locations of the defects and 5 C is the covariant 
version of the Dirac delta function given by 

S c (a,a')= lim A eX p (-~<P(a, a')) = 5 2 (a - a')/^/g(a) 

A^oo 2,71 \ I J 

(d(<r,a') is the geodesic distance between a and a'). We note that the Gaussian 
curvature K = K(cr) plays the role of a background charge density. Because of 
this, positive defects tend to concentrate in regions with positive curvature, whereas 
negative prefer saddle shaped regions with negative curvature. Charge neutrality and 

\ The Greens's function of the negative of the Laplace-Beltrami operator A is the "electrostatic" 
potential at the point a of a unit charge located at point a' and of a negative unit charge which 
is uniformly distributed over the surface to ensure charge neutrality. It is given by — AG(a, a 1 ) = 
5 c (o, a') — l/A, where S c is the covariant delta function and A is the area of the surface. 
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the Gauss-Bonnet theorem |D| 2ir ^\ Qi — J dA K = yield the topological constraint 
E,?i = 2(1-7)- 

Above a certain temperature it is expected that the low-temperature phase 
with a few tightly bound defects is destroyed through unbinding of defect pairs, 
analogous to the Kosterlitz-Thouless transition in flat space. The high-temperature 
phase has a finite density of thermally excited, unbound defects. The interaction 
between the defects is screened, with a screening length of the order of the typical 
distance of the defects (Debye-Hiickel length). Above the transition temperature, we 



make a Gaussian approximation of the Coulomb gas model 2/2 with a continuous 
defect density p 

Y = ^rJ dA J dA ' ~ K) ° G(a ' a>) (2 ^° ~ K) °' + hj dA p2 ' (2 - 3) 

where x is the fugacity of the charges. By setting SH/Sp = 0, we obtain for the mean 
charge density 

2n P = ,_\^^ = K + + 0(x->). (2.4) 

Although this approximation accurately represents the Coulomb gas at high 
temperatures, the use of the Coulomb gas at all to describe the high-temperature 
phase is rather suspect. Nevertheless, we show that the Coulomb gas model yields a 



density of defects which agrees remarkably well with the density obtained from (2.1) 
in the high-temperature phase on an arbitrary curved surface. 

3. Charge density in the high-temperature Gaussian approximation 

In the remainder of the paper we present the calculation of the defect density p from 
( [2.l| ) in the disordered state, where the quartic term J dAiuiU 1 ) 2 is irrelevant and can 
be neglected. We expect that similar to the situation in curved space-time Q , a term 
proportional to J dA K UiU 1 is generated under renormalization, where K — K(a) 
is the Gaussian curvature. Since K has the dimension of 1/lcngth 2 this term is as 
relevant as the gradient term. Thus, in the high temperature phase the vector field is 
distributed according to the Gaussian weight 



P[u] oc exp f-i J d 2 o^fg(DiU D l u j 



+ TU l {a)ui{a) + J] Ku l Ui) , (3.1) 



where r is now the mass of the vector field and r\ is the coupling of K to UiU 1 . In 
addition, the distribution for U{ has to be equipped with a covariant cutoff procedure, 
such as the heat kernel regularization [|l9). Because the model depends only on the 
intrinsic geometry of the manifold, no extrinsic couplings (such as a term proportional 
to C 2 UiU z , where C is the mean curvature of the surface) can be generated under 
renormalization. 

The zeros of the field Ui are characterized by the index q = signdet(Di7ij) = ±1. 
The index (charge) describes the local topology of a flow Ui near a zero u, (ct) = 0. 
The corresponding mean charge density is given by 

p(0 = (Z)9i*c(^(r i )) ) (3.2) 
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where the defects are labeled by the index i, located at coordinates <7j, and have 
charges q^. The expectation value is taken with respect to the probability distribution 
of ©. 

Transforming from the variable a to the variable u via the Jacobian, we obtain 

p = (det(Aw J ' O)) 5 C (u(tr))) , (3.3) 

which can be seen easily using a locally Euclidean coordinate system and linearizing 
the vector field around the zero Uj(x\,X2) — x k a k f. 

d 2 x det {di(x k a kj )) S 2 (x k a kj ) 

= J d 2 x | det(a,j)| sign det (ajj) 5 2 (x k a k j) — sign det (ay ) = ±1. (3.4) 

To calculate the expectation value (at point a) p{cr) one needs the joint distribution 
of Ui(a) and DiUj(a) which can be determined since Ui(a) and DiUj(a) are a 
set of (multicomponent) Gaussian random variables with correlations (ui(cr)uj(cr)) , 
(u k (a)DiUj (a)} , (DiUj (a)D k u t {a)}. 

3.1. Density of defects in fiat space 

Before calculating results in curved space, we review Halperin's calculation for the 
density of zeros of a Gaussian two-component order parameter u in two dimensional 
flat space ||. Equation (3.1)) becomes 

P[u(r)} ocexpj-i J d 2 r [{d^) 2 + tu 2 ] J (3.5) 

and (3^) becomes 

p(r) = (6 2 [u(r)}det(d l u J )). (3.6) 

This expectation value is completely determined by the probability distribution 
P(£i, ctij), where £j = Ui(r) and = diUj(r) via the formula 

p(r) = f d 4 otijP (o, a>ij) det ctij. (3.7) 



Since (3.5) is Gaussian, this probability distribution is just given by 



1 1 f 1 



P&W) = ^——^^--xtMZ x,| , (3.8) 

where a; is a six-component vector given by x — (^i, £2, cti2> &2i, 0*22), and My is 
the matrix of correlations Af,j = (xiXj). Plugging (3.8) into (3.7) gives 

1/2 



P= 7T- 



det Mij 
2vr \ det My 



[m u ~ M 23 ) , (3.9) 



where M is the matrix of correlations My = (t/jj/j), and y is a four-component vector 
given by y = (an, a i2 , a 2 i, a 22 ). 

The expectation values necessary to evaluate (3.9) can be readily determined 
from (|3.5|). The result is that p(r) = 0, as expected by symmetry: the system is 
uniform and charge-neutral. In order to obtain a nontrivial result to compare with 
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the Coulomb gas model, we must calculate the correlation function of the density of 
charges, defined by 

C(r) = (p(r)p(o)) = (S 2 [u(r)} det S 2 [u(o)} det [$^(0)]) . (3.10) 

This can be evaluated by similar methods, and the results match up very well with 
the Coulomb gas model Q. 

3. 2. Density of defects in curved space 

In the case of a general closed membrane, p will already be nontrivial for two reasons. 
First, the system is not charge-neutral, but rather the total charge must be equal 
to two minus the number of handles on the surface (e.g. two for a sphere, zero for 
a torus, etc.) Second, unless the surface has a high degree of symmetry, the charge 
will not distribute itself uniformly. Rather, the charge density will depend upon the 
local curvature of the surface. Thus, for the case we consider in the remainder of the 
paper, we can get a meaningful comparison with the Coulomb gas theory solely from 
calculating the charge density p, rather than needing to calculate the more complicated 
correlation functions. 



The method used is conceptually the same as in section 3.1, but more technically 
complicated due to the curvature of the space. Therefore, we present it in Appendix 
|A|. The analogous result to (3.9) for a general Gaussian, 0(2) invariant distribution 



for the vector field it, is 



2np K = e*V ! .Dfc " J '.' , (3.11) 
V {u m u m ) 



Since the right hand side of (3.11) is a total divergence, we find, after integration 
over the surface the aforementioned topological constraint for the total charge of the 
defects 

2tt J dA p = J dAK = 4tt(1 - 7) (3.12) 
which agrees with the Gauss-Bonnet theorem, as 1 — 7 is the genus of the surface. 



To derive p from ( 3.11 ), we need to calculate (uiUj) and {(DiUj)uk). This can 
be done in an expansion with respect to the interaction range 1 /r using the Gaussian 
weight ( |3.l[ ) . It is convenient to decompose the vector field u, into a sum of a gradient 
and a curl m = di<j) + ti 3 dj\. This representation is only valid for (deformed) spheres. 
For other topologies, modes exist which cannot be written as sum of a gradient and 
a curl. For example, in a torus, a vector field that represents a flow along one of the 
perimeters cannot be decomposed in this way. 

A particularly simple case is given for 77 = 1 because the potentials <j) and x 
decouple. The special role of the 77 = 1 case can also be understood within the gauge- 



field representation, as shown in Appendix D. For this case, we derive the density of 



defects in a high-temperature expansion. For high temperatures, the screening length 
is small compared to the radius of curvature: the surface appears to be almost flat. 
Upon increasing screening length, more details of the geometry become relevant. We 



present the details of this high-temperature expansion in Appendix B, obtaining as 
our main result the average defect density p: 

where Z 1 = 2tt/ log(l/(o 2 r)), and a is a short distance cutoff. 
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To lowest order in the correlation length t -1 / 2 , this is equivalent to the Debye- 
Hiickel approximation ( |2.4D provided one identifies irK^x — 3Zt. For larger 
correlation lengths, however, deviations show up. The term oc r _1 will be independent 
of the coupling r/ for dimensional reasons. The next orders, however, depend on rj. We 
conjecture, that the expansion (3.13) remains valid for arbitrary genus of the surface. 
Treating general genus and 77, however, requires the calculation of moments of the 
vector field ui directly, which is much more complicated and beyond the scope of this 
work. 

It will be difficult to observe the defect density ( 3.13 ) experimentally, since p, 
which is density of positive defects minus the densi ty of negative defects is of the order 
of p ~ 1/A due to the topological constraint ( 3.12 ). On the other hand the density of 
positive defects plus the density of negative defects is of the order of l/£ 2 ~ t, where £ 
is the correlation length, which is small well above the transition temperature. In the 
high-temperature phase, therefore, we have to measure a tiny density difference in the 
presence of a large background density. Closer to the transition region the background 
density becomes smaller and there might well be a chance to resolve the defects in 
thin films using polarized light. It would be certainly interesting to see wether our 
expansion ( |3.13| ) or the result (2.4) obtained from the Coulomb gas model allow for 
the better fit to the experimental data. 



4. Conclusion 



We have derived the average topological charge density of vector fields with a Gaussian 
distribution on a curved surface. We found that for high temperatures, the zeros 
behave like (screened) charges in the presence of a background charge density equal 
to the Gaussian curvature. We demonstrated the validity of the Debye-Hiickel 
approximation of the Coulomb gas model, which, as discussed in section [l], is not 
obvious, since the Coulomb gas model originates from a low-temperature model of the 
orientational order and we are attempting to apply it in a high-temperature regime. 
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Appendix A. Calculation of charge density 

We will calculate the density of defects p = (de^DiU 1 * (a)) 6 C (it(<r))). We begin by 



exact analogy with the calculation for flat space outlined in section 3.1, by noting 
that this expectation value is completely determined by the probability distribution 
P(vi, Aij), where Vi = Ui(a) and Aij = DiUj(a) via the formula 

p{a) = / d ; .l,,/'u>..l,.,:<lrt .1,'. (A.l) 
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(A.2) 



Since (3.1) is Gaussian, this probability distribution is just given by 

where a;, is a six-component vector given by (i>i, «2, An, A\2, A21, A22), and Mjj is the 
matrix of correlations My = (xiXj). 
We can reexpress this as 

P(«i, ^) = ^ y d 4 i^ J d 2 v l exp j-^PM^ + is^ij , (A.3) 

where J* is a six-component vector g iven by (v 1 , v 2 , A 11 , A , A 21 , A 22 ). Performing 
the integral over the x % returns us to ( A.2 ). Explicitly we have for the defect density 

P = J d i A^d 4 A lJ d 2 Vd 2 v l exp (-^A kl ((A«i)(Pk«j))) 

x exp ^-A i:i v k ((DiUj)u k ) - ^v l v ] (uiUj) + iA tJ A tj + itfv^j 

x^e ik f? l A i3 A kl 5 c {v). (A.4) 
The integration over the Vi yields a factor of y/g since 

<*c (v(o-)) = firn A exp (~g ij (<^>iV^j = y/g(a) 6 2 {v). 

The resulting expression can be simplified by using the 0(2)-invariance of the field Ui 
by noting that is the only rank two tensor invariant under 0(2), and therefore 

(UiUj) = (u m u m ) . (A. 5) 

After integration over the v l we obtain 

P 



n (u m u m ) (2tt) 4 J V 2 / 

t I en \(n \\ ((AMj> m ) {{D k uj)u m ) 
Tijki = {{D t Uj){D k ui)) - 2 — . (A. 7) 

(u n u ) 



2 

where 



Integrating over the A l i in flA.q ), we see that the Ai 3 are simply Gaussian variables 
with correlations given by 

(AijAki) =T ijkl . (A.8) 

Therefore, JA^ ) yields 



2tt (u m n m ) 

To further simplify this equation, we derive 

jkjl n ( {{D l Uj)ui) \ ik . t {{D k D t u 3 ) Ul ) ikJl ({D z u 3 )(D kUl )) 
e e J L) k I — ; — — I = e e J ; — he e 



(A.9) 



u m u" 1 ) 1 {u m u ) [u m u" 



2e ik £ ji ((A"j)m) ((£>fcU m )M m ) 



u n u"-) 2 



■7. _'l 
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using the fact that in two dimensions, the Riemann curvature tensor is Rijki — tij^kiK, 
where K is the Gaussian curvature), and also that due to the 0(2)-invariance of itj, 

((A%H) = \g 3k ((D lUm )u m ) + \e ]k e mn ((D t u m )u n ) . (A.ll) 
We therefore arrive at an expression for the mean defect density 

2, P -K = e^D k (^m, (A.12) 



(u m u m ) 

valid for any Gaussian, 0(2)-invariant distribution for the vector fields Ui. 
Appendix B. high-temperature expansion 

In this Appendix, we derive the densi ty o f defects in a high-temperature expansion for 
the case 77 = 1, as disucssed in section 3^. In this case, since D l Didj(j>—Kdjcj) — djAcj), 
the eigenfunctions of the operator —DiD 1 + K + r (acting on vector fields) 

{-D l Di +K + t) Uj, a = (A Q + t) u,, a (B.l) 

can be written as = di4> a and u^l — e^dj^a, where 4> a is a normalized 

eigenfunction of the Laplace-Beltrami operator —A(f> a — —g 13 Didj(f> a — X a <p a . 
Together with the normalization of the u a 



dA Uj^ a u 0,a — J dA g v di(f> a dj(j) a = X a , 
we obtain the propagator for Ui 

(ui(a) Uj (a')) = (d t <P(a)d 3 <P(a')) + e; fe (a)e/ (a') (d k ^(a)di<P(a')) (B.2) 
in terms of the scalar propagator 

{<K<r)4>{o>)) = E Tit^rr M^M**), (b.3) 

^— ' A Q (A Q +r) 

where the zero mode <j) =const is omitted and a is an exponential cutoff length that 
arises from the heat kernel regularization. After some algebra we find 

2wp - K = e e D k { j = -D ^ j , (B.4) 

with ((A(f>)d t <j)) = (l/2)di((A4>)(j)) from ([rJ|). Both - {<t>A<j>) and (d m <f>d m 4>) are 
logarithmically divergent for small cutoff lengths 

{ - m {tZ) ] } = h log (i) + fimte parts ' 

where only the finite parts depend on the position on the manifold. Therefore the 
numerator of (B.4) is finite, whereas in the limit of small cutoff lengths the divergent 
denominator can be replaced by its most divergent (spatially constant) part. Defining 
Z- 1 = 2tt/ log(l/ (a 2 r)), we obtain 

2np-K = Z^A (((f)A<f>) ~ (d m <pd m cj))) . (B.5) 

Using (d m (/>d m (j)) = (1/2) A (<f) 2 ) - (0Ac/>) and t (0 2 ) = (<j>A<f>) + (4>(t - A)0>, we 
obtain 

2np - K = Z- 1 U (2 - ±a\ (0A0) - ^A 2 (0(r - A)0)) . (B.6) 
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The moment 

{4>{t - A)(p) = > 4> a (ct) 

K 

is the Green's function of the Laplace-Beltrami operator at coinciding points, which 
can be obtained from conformal field theory. With the definition of the massless 
propagator and its short distance behavior G(er, a') ~ r(er) — \og{d{a, a'))/(2n), we find 
for two conformal equivalent metrices gij = (g^ that — Ar + 2/A — K/(2n) = (-AT + 
2/ A — K/(2it))/( and consequently for a spherical topology — Ar = K/(2n) — 2/ A. 
A is the area of the surface and d the geodesic distance. We have A 2 (0(r — A)<j>) = 
—AKJ (2tt). The moment (4>A(f>) has to be calculated in an asymptotic 1/r expansion. 
With the help of the references (19, 20, |T) we find for the finite part of (4>A(f>) 



K K 2 + AK ., 



- (<^) = — + 2 + 0(r- d ). (B.7) 

127TT o(J7rr z 

Finally, we obtain the average defect density p 

— ^ + «-3^+°<^ <■»> 

Appendix C. The gauge- field representation 

Orientational order is frequently represented as a gauge-field theory J^, || . The vector 
field Ui is represented in a local orthogonal base (reference frame) Vi by a complex 
function ij) through m — u,: Re("0) + e^Vj Imhp), where ViV 1 = 1 and e is the 



antisymmetric unit tensor. Plugging this into (2.1), we obtain 
H 



T =J d 2 aV5 {9 l] (di^ + W) (drf - iQ^) + t\iP\ a + C |V| 4 ) (C.l) 

with the vector potential fi^ = e^VjDiVk, resulting from the fact that the reference 
frame Vi is changing from point to point. In this representation the Gaussian curvature 
K plays the role of a perpendicular magnetic field e^DjOj = K. Since any unit vector 
field Vi must have two points on a sphere where it is singular, will be singular at 



two points as well, even if the underlying surface is not. We therefore use (2.1) rather 



than (C.l) in this paper as a base for calculation. The disadvantage of (2.1) is that 
it is more difficult to generalize to an n-fold symmetry, as is done by Park, et al 
and Evans ||. 

Appendix D. The r\ = 1 case 



Equation (C.l)) is covariant and gauge invariant, i.e. invariant against changes of the 



local frame Uj. For convenience, we choose a conformally flat coordinate system with 
a metric tensor = £(x, y)Sij, where the coordinates are a 1 = x, a 2 = y and £(x, y) 
encodes the (intrinsic) geometry of the surface. A proof that such a coordinate system 
exists, as well as a review of its properties, can be found in (Note, however, that 
for closed surfaces other than tori, Q will have singular points.) Furthermore, we 
choose a particular reference frame v x = yf, v y = 0. Then ^fgg 1 ^ = Sij, il x — d v u> , 



il y = -d x cu, K = (1/C)(<9 2 + d 2 )uj where uj = -(1/2) logC The Gaussian weight fl3.l| ) 



becomes P[ip] = cxp(— 1/2 J d 2 x tfj*Hip) with the Hamilton operator 

H = -(d x - idyLo) 2 - {d y + id x uj) 2 + r/V 2 w + Qt. (D.l) 
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For 77 = 1 and a vanishing mass r = we can express the Hamilton operator in terms 
of the square of a Dirac-type operator [a x VjZ are the Pauli-matrices) 

H = — [(d x - idyUj) a x + (d y + id x uj) cr y ] 2 

= -(d x - idyUj) 2 - {d y + id x uj) 2 + a z V 2 uj (D.2) 

in the a z = +1 sector. The latter operator is the Hamilonian of the 2D-Pauli equation 
for spin— 1/2 particles with the (dimensionless) magnetic moment g = 2 (electrons!) 
in a magnetic field B z = V 2 uj. Thus the 77 = 1 case that we focus on is closely related 
to both the Dirac equation and the Pauli equation in two dimensions. The Dirac 
equation and its discretized counterpart are commonly used to describe the properties 
of electrons confined to a plane in a quenched, perpendicular magnetic field |p2[ . 
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Figure 1. A surfactant molecule tilted away from the normal. 
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Figure 2. The planar vector field (u x ,u y ) = (x 2 — 1,2/) with a negative zero 
(left) and a positive zero (right). 



